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ABSTRACT 

We have studied the dynamics and masses of a sample of ten nearby luminous and ultraluminous 
infrared galaxies (LIRGS and ULIRGs), using 2.3 /im ^■^CO absorption line spectroscopy and near- 
infrared H- and i^T^-band imaging. By combining velocity dispersions derived from the spectroscopy, 
disk scale-lengths obtained from the imaging, and a set of likely model density profiles, we calculate 
dynamical masses for each LIRG. For the majority of the sample, it is difficult to reconcile our mass 
estimates with the large amounts of gas derived from millimeter observations and from a standard 
conversion between ^^CO emission and II2 mass. 

Our results imply that LIRGs do not have huge amounts of molecular gas (10^° — 10^^ M©) at 
their centers, and support previous indications that the standard conversion of ^^CO to II2 probably 
overestimates the gas masses and cannot be used in these environments. This in turn suggests much 
more modest levels of extinction in the near-infrared for LIRGs than previously predicted (Ay ^ 10- 
20 versus Ay ~ 100-1000). The lower gas mass estimates indicated by our observations imply that 
the star formation efhciency in these systems is very high and is triggered by cloud-cloud collisions, 
shocks, and winds rather than by gravitational instabilities in circumnuclear gas disks. 
Subject headings: galaxies: kinematics and dynamics, galaxies: starburst, galaxies: stellar content 



1. INTRODUCTION 

Luminous and ultraluminous infrared galaxies (LIRGs 
and ULIRGs) have extremely high bolometric luminosi- 
ties and have spectral energy distributions that are 
dominated by infrared light (8-1000 /^m luminosities > 
10" Lq for LIRGs and > lO^^ ^q for ULIRGs). Their 
morphologies are generally disturbed, sometimes with 
tidal tails or two discernable nuclei, serving as strong 
evidence that they are the products of mergers and in- 
teractions (e.g., Farrah et al. 2001; Bushouse et al. 2002; 
Veilleux et al. 2002). Their infrared power is thought to 
be from intense star formation triggered by the merging 
process, supplemented by the outputs of active galactic 
nuclei (AGNs). It has been suggested that such systems 
go through not only a luminous starburst phase, but later 
evolve into QSOs (Sanders et al. 1988). 

LIRGs are often the result of the collision of two gas- 
rich disk galaxies and have been suggested to be the pre- 
cursors to elliptical galaxies (Toomre & Toomre 1972; 
Schweizer 1986; Hernquist 1992, 1993). This hypothe- 
sis can be tested using the fundamental plane, a three- 
parameter relation among effective radius, velocity dis- 
persion and mean surface brightness for early-type galax- 
ies (e.g., Djorgovski & Davis 1987; Dressier et al. 1987). 
It has been shown (Genzel et al. 2001; Tacconi et al. 
2002) that most LIRGs are remarkably close to the fun- 
damental plane of dynamically hot galaxies, particularly 
on the less "evolution-sensitive" effective radius vs. ve- 
locity dispersion projection. Only LIRGs with the small- 
est effective radii are offset from the plane. The place- 
ment of LIRGs within the fundamental plane matches 
the locations of intermediate-mass, disky (rotating) el- 
lipticals. The closeness of LIRGs to ellipticals on the 

^ Observations reported here were obtained at the MMT Obser- 
vatory, a joint facility of the University of Arizona and the Smith- 
sonian Institution. 



fundamental plane is somewhat surprising, but is usu- 
ally interpreted as a "conspiracy" of stellar evolution and 
extinction - the surface brightening of the new stellar 
population is obscured by heavy dust extinction to keep 
the LIRGs on the plane. However, the dynamical sim- 
ilarities of the LIRGs and intermediate-mass ellipticals 
are robust and do not depend on the photometric prop- 
erties. Therefore, LIRGs are likely an important link 
to understanding the formation and evolution of early- 
type galaxies from gas-rich spirals. This link is now be- 
ing strengthened by studies of the mid- and far-infrared 
properties of LIRGs and ULIRGs at high redshifts (Char- 
mandaris et al. 2004; Egami et al. 2004; Le Floc'h et al. 
2004) , by ground-based integral field spectroscopic stud- 
ies investigating the dynamics of low redshift systems 
and their link to deriving dynamical masses for LIRGs 
at higher redshifts (Colina et al. 2005), and by spectro- 
scopic expansions on Genzel et al. (2001) and Tacconi et 
al. (2002) studying dynamical evolution (Dasyra et al. 
2006). 

In the conventional view of LIRGs, one by-product of 
the merger of two gas-rich disks is that huge masses of 
molecular gas will settle into their nuclei. Interferomet- 
ric observations of LIRGs at millimeter wavelengths in- 
dicate large concentrations of molecular gas at their cen- 
ters. Calculations to derive gas masses using maps of the 
^^CO J= 1^0 transition sometimes adopt a standard 
conversion factor between the ^^CO J = 1 ^ luminosity 
and the molecular hydrogen (H2) mass. The conversion 
was originally derived from observations of ^^CO to II2 
gas in giant molecular clouds in our own galaxy (Scoville 
& Good 1989) and takes the form 
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where Mg is the gas mass (Af(Il2 + He)), Fco is the 
spatial- and velocity-integrated CO flux in Jy km s^^. 



Dmjic is the distance to the galaxy in Mpc, and z is the 
redshift. The standard conversion factor is equivalent to 
N{H2)/Ico = 2.24 X 10^0 cm-^ (K km s'^)'^ . This 
conversion produces molecular gas masses in the range 
10^° - lO"Af0 (Bryant & Scoville 1999), indicating ex- 
tinctions as large as Ay ^ 100 — 1000 in the optical 
(see, e.g., Narayanan et al. 2005 for recent alternative 
conversion factors). If true, this would mean that even 
observations taken in the near-infrared would barely pen- 
etrate the surface of most LIRGs and that much of their 
star formation activity would remain hidden from view. 

However, there are indications that the standard con- 
version between CO flux and H2 mass may not be ap- 
plicable in the environment in LIRGs. Both theory and 
observations show that the state of gas in the centers of 
LIRGs is quite different from that in local giant molec- 
ular clouds. Maloney & Black (1988) argued that the 
high densities and temperatures of the molecular gas in 
dense starbursts would result in the standard ratio sig- 
nificantly overestimating the total mass. They stated 
that, while the conversion does seem to be constant over 
most of the Milky Way, the underlying reason for this 
is not well understood, and should not be assumed for 
other galaxies. Additionally, a CO-to-H2 ratio that varies 
across the Milky Way presents a solution to the apparent 
observed discrepancy between the radial gradient in the 
diffuse Galactic 7-ray and the distribution of supernova 
remnants (Strong et al. 2004). There is also evidence for 
variances in the standard conversion between galaxies, 
where arguments for a strong dependence on metallicity 
as traced by the oxygen abundance [0]/[II] have been 
presented (Israel 1997; 2000; Barone et al. 2000; Boselli 
et al. 2002). 

Solomon et al. (1997) and Downes & Solomon (1998) 
concluded that the standard factor may be high by a 
factor of five based on a comparison of dynamical and gas 
masses derived from far-infrared observations. Lisenfeld, 
Isaak, & Hills (2000) also suggested, using submilhmeter 
measurements of cold dust emission, that the molecular 
gas mass deduced from CO observations may be 2-3 times 
higher than would be deduced from thermal emission. 

Other works have challenged the molecular gas masses 
calculated from the standard conversion by comparing 
those values to measurements of the total dynamical 
masses. Shier, Rieke, & Rieke (1994) obtained dynami- 
cal masses for three LIRGs based on velocity dispersions 
from near-infrared spectra and scale-lengths derived from 
near-infrared light profiles. They found that in all three 
galaxies the total dynamical masses of the cores were 
less than the claimed H2 masses. They concluded that, 
despite large errorbars on the dynamical masses, the 
standard conversion factor must be wrong for at least 
some luminous infrared galaxies. Alonso-Herrero et al. 
(2001) calculated a dynamical mass for NGC 1614 de- 
rived from Br 7 spectroscopy (Puxley & Brand 1999) and 
their knowledge of the morphology and size of a ring of 
star formation revealed in _ffS'T/NICMOS Pa a images. 
They found that the total dynamical mass of the star- 
burst region is nearly four times smaller than the mass of 
molecular gas estimated from the standard ratio of ^^CO 
(1 ^ 0) to H2. 

The difference in estimated masses is important for un- 
derstanding the way in which stars are formed in these 
galaxies. For instance, if the large surface densities of 



gas in the circumnuclear disks predicted by the standard 
ratio are correct, then star formation occurs through in- 
stabilities in the disk (Taniguchi & Ohyama 1998). How- 
ever, if lower mass estimates are valid, then star for- 
mation through cloud-cloud collision, shocks, and winds 
(Scoville, Sanders & Clemens 1986) is predicted. Ad- 
ditionally, using the standard conversion suggests a low 
star formation efficiency, whereas lower mass estimates 
yield high star formation efficiency. 

To understand more fully star formation and extinction 
in LIRGs, we determine dynamical masses for a sample 
of ten galaxies. We describe the observations and data 
reduction in ^ discuss the analysis and modeling in '^ 
compare the molecular masses to dynamical masses for 
the sample in ^ and summarize the results in |S| 

2. OBSERVATIONS AND DATA REDUCTION 

Ten luminous and ultraluminous infrared galaxies were 
chosen for this study. These galaxies have K-h&Tid mag- 
nitudes in the range ~ 12-13.5. All of the LIRGs are 
nearby {v < 17,000km s""'^), and many have millimeter 
CO emission measurements (e.g., Sanders, Scoville, & 
Soifer 1991) as well as NICMOS-2 images of their inner 
regions (Scoville et al. 2000; Alonso-Herrero et al. 2000; 
Alonso-Herrero et al. 2001). Thus, we are able to choose 
galaxies with very compact nuclear CO associated with 
high molecular gas mass and with simple nuclear struc- 
ture (usually an r^'^ or exponential profile). 

Near-infrared images {H and Kg) of the majority of 
galaxies in the LIRG sample were obtained with the 
wide-field PISCES camera (McCarthy et al. 2001) at the 
2.3 m Bok and 6.5 m MMT telescopes. Table 1 gives the 
dates and exposures times for each galaxy in both bands. 
All images were corrected for quadrant cross-talk effects 
known to be present in 1024 x 1024 HAWAII arrays, 
using a custom C program written by Roelof de Jong. 
The images were then dark-subtracted with a combina- 
tion of 11 dark frames of the appropriate exposure time. 
Each frame was divided by a normalized flat-field cre- 
ated by a median combination and sigma-clip rejection 
of the dithered science frames. Bad pixels were masked 
out using a script that substitutes known hot or dead 
pixels with the value of the mode of the image. The 
PISCES field-of-view is an inscribed circle, so the outer 
edges of the images were also masked. Frames were cor- 
rected for geometric distortion, although PISCES' dis- 
tortion is known to be less than a 3% effect at the field 
edge and so does not significantly contribute to poor im- 
age quality or errors in photometry. Images were then 
aligned and stacked using standard IRAF tasks. Seeing 
was ^ 1'.'2-1'.'5 for the images taken at the 2.3 m and ^ 0'.'6 
for images taken at the MMT. For two of the galaxies, 
NGC 1614 and IC883, we rely solely on the NICMOS-2 
observations and do not have wide-field images. i/-band 
images for the entire sample are shown in Figure 1. 

The morphologies of the galaxies vary from late-type 
spirals to irregular mergers and interactions (see Ta- 
ble 1). Tidal tails and other features indicating recent 
merger activity can be seen in several of the LIRGs (e.g., 
NGC 1614). 

Near-infrared spectra were taken at the MMT on 2000 
December 5 and 2001 March 2-3 with the FSpec spectro- 
graph (WiUiams et al. 1993). The sht width was ~0'.'7, 
corresponding to two pixels. The seeing- limited images 



were smaller than the slit width for all the observations. 
The 600 lines mm~^ grating was used at r^ 2.3 /im to ob- 
serve the first overtone CO absorption bandheads, where 
one pixel corresponds to ~ 3.3 A. Table 1 shows the ex- 
posure times for each galaxy. The terrestrial absorption 
spectrum was sampled by observing a star of type A V 
- G V between each separate galaxy observation. Cross- 
correlation template stars of K5 III - M2 III were also 
observed. 

All spectra were reduced as described by Engelbracht 
et al. (1998). Frames were dark-subtracted, trimmed, 
masked for bad pixels, sky-subtracted, and flat-fielded 
using IRAF tasks called by a custom script. Multiple 
frames of objects were aligned and combined, and one- 
dimensional spectra of objects, standard stars, and sky 
frames were extracted. The one-dimensional galaxy and 
standard star spectra were wavelength calibrated using 
the wavelengths of OH airglow lines. Each galaxy spec- 
trum was then divided by a template standard star spec- 
trum of type K5 III - M2 III. The continuum in each star 
and galaxy spectrum was fitted with a low-order polyno- 
mial that was then subtracted. The spectra are shown 
in Figure 2. 

The first overtone CO absorption features are one- 
sided, but sharp and strong and hence suitable for studies 
of dynamics. The measured bandheads have rest wave- 
lengths of 2.2935 /xm for ^^CO J= 2 -* and 2.3227 ^^m 
for 12C0 J= 3 ^ 1 (Kleinmann & Hall 1986). The first 
CO overtones are filled in by emission by hot dust for Mrk 
231. Therefore, for this galaxy we measured the second 
CO overtones at 1.67 /im, using similar procedures for 
both the observations and the reductions as for the rest 
of the sample. 

3. ANALYSIS AND RESULTS 

3.1. Velocity Dispersions 

To determine dynamical masses for the LIRGs, the 
velocity dispersions must be calculated from the near- 
infrared spectra using template stars with similar CO 
absorption bands. There are three main techniques that 
can be used: the Fourier quotient method (Sargent et 
al. 1977), the cross-correlation method (Simkin 1974; 
Tonry & Davis 1979), and direct fitting (Franx, lUing- 
worth, & Heckman 1989; Rix & White 1992). In the 
Fourier quotient method, the galaxy is assumed to be the 
convolution of a single stellar template with a broaden- 
ing function, B, which can be retrieved using an inverse 
Fourier transform. The disadvantage of this method is 
that the errors in the quotient are correlated, so that er- 
ror analysis becomes difficult. Also, because the broaden- 
ing function is fitted in Fourier space, absorption features 
from all parts of the spectrum interact with each other, 
making B very sensitive to template mismatches. In the 
cross-correlation method, the galaxy and template are 
correlated directly, and the location of the highest cor- 
relation peak is identified with the mean velocity. The 
dispersion is estimated by subtracting the width of the 
template autocorrelation peak in quadrature from the 
width of the cross-correlation peak. The drawback of this 
method is that the spectra must be padded with extra ze- 
roes, possibly with some trimming of the spectral edges, 
resulting in lost information. The direct fitting method 
of Rix & White (1992) assumes that the galaxy differs 
from the template by the shape of the continuum and 



the kinematic broadening. The broadening is thought of 
as a superposition of multiple templates shifted in ve- 
locity space. The disadvantage for this method is that 
numerous templates of differing stellar types are almost 
certainly needed for a good fit to the galaxy spectrum, 
taking up more observing time. Also, this type of fitting 
program may get trapped in local minima if the initial 
guesses for radial velocity and velocity dispersion are not 
reasonably accurate. 

Of these choices, we used the cross-correlation method 
on our spectra. Although the spectra have to be care- 
fully continuum-subtracted and padded with zeroes, we 
have trimmed very little information from either edge, 
and have tested that losing a few Angstroms of spectrum 
does not change the velocity dispersion results within the 
calculated error bars. 

Velocity dispersions of the galaxies were calculated by 
comparing the width of their CO absorption bandheads 
to those of the template standard stars. We first gen- 
erated a series of broadened template star spectra with 
Gaussian dispersions corresponding to a range in veloc- 
ity of 20-200 km s^^. The broadened spectra were then 
cross-correlated with the original template star spec- 
trum, and the width of the cross-correlation peak as a 
function of the Gaussian broadening was measured. Each 
galaxy spectrum was then cross-correlated with a tem- 
plate star. The width of the cross-correlation peak be- 
tween the galaxy and template star revealed the velocity 
broadening, based on the broadened stellar template cal- 
ibration. Velocity dispersions and uncertainties for each 
LIRG are listed in Table 2. Uncertainties were estimated 
by looking at the range of dispersion values obtained from 
comparing multiple template stars with each galaxy. 

Mkn231 was observed at 1.67 /xm to calculate the ve- 
locity dispersion from the CO bandheads there. How- 
ever, due to difficulty in subtracting the continuum at 
this wavelength range, we used the Tacconi et al. (2002) 
value listed in Table 2 for the remainder of this work (see 
Shier 1995 for a discussion of continuum subtraction). 

Table 2 shows a comparison of our dispersion mea- 
surements to others. Three of our galaxies, NGC 1614, 
NGC 2623, and IC 694/NGC 3690 ^, overlap with a study 
by Shier, Rieke, & Rieke (1996). We find that the val- 
ues for the velocity dispersions are higher for the first 
two than those found by Shier et al. (1996), while the 
value for NGC 3690 is in agreement. One possibility for 
the differences is the differing areas of the galaxies sam- 
pled by us at the MMT (slit width - 0'.'7) and by Shier 
et al. at the Bok 2.3 m (slit width - 2'.'4). Where the 
CO band strength strongly peaks on the galaxy nucleus, 
our measurements are strongly weighted by the center of 
the LIRG and may overestimate the dispersion. Alonso- 
Herrero et al. (2001) show that the CO photometric 
index in NGC 1614 is extremely large in the inner 0'.'4 
in comparison to the Pa a emission. Thus, an aperture 
effect is a likely explanation for the discrepancy on this 
galaxy. For NGC 2623, it is also possible that the differ- 
ence in slit width is accounts for the (smaller) discrep- 
ancy in dispersion measurements, although we do not 
have the benefit of high resolution NICMOS CO imag- 

^ There has been some confusion in the literature regarding 
which galaxy is actually IC 694; we refer readers to Hibbard & 
Yun (1999) for a discussion. 



ing to provide supporting evidence. Within the uncer- 
tainities, our dispersion values for IRAS 10565+2448 and 
Mkn273N agree with the values derived by Dasyra et al. 
(2006) and Tacconi et al. (2002), respectively 

3.2. Light Profiles 

To obtain dynamical masses for the sample, bulge 
and/or disk scale-lengths must be fitted to the large- 
scale near-infrared images. Two-dimensional models of 
the galaxies were created using a de Vaucouleurs, r^^^ 
bulge and exponential disk fitting routine (Rix & Zarit- 
sky 1995). The routine rebins the images onto a cylin- 
drical grid, given a specified minimum and maximum 
radius and the coordinates of the galaxy center. Fed an 
initial guess for a set of bulge and disk parameters, the 
program begins a "biased random walk" within a pa- 
rameter box until it finds the combination of parameters 
that minimize chi-square. Typically, a few pixels at the 
very centers of the galaxies are not included in the fit, 
because the fitting routine does not account for the see- 
ing. Although both H and Ks-ha,nd images were fitted 
where available, we report the results for iJ-band, where 
we have images for the whole sample. Exponential disk 
and r^'^ effective bulge _ff-band radii with errorbars are 
listed in Table 2. 

3.3. Moss Models 

As LIRGs are the products of mergers or are mergers in 
progress, it is not clear whether the inner regions of the 
galaxies should be expected to be hot spherical systems 
or cold disk systems. Radial profile fits to a large sample 
of LIRGs (e.g., Scoville et al. 2000) show that the best 
fitting models can be r^'^, exponential, or neither. If the 
LIRGs are ellipticals in formation, then we might expect 
their centers to be kinematically hot spheroids. However, 
young stars with strong GO lines are formed in the gas 
that falls to the center of the galaxy and may have settled 
into a ring or disk (Barnes & Hernquist 1991). Also, 
a remnant disk consisting of the old stellar population 
of an original merging spiral galaxy may still be intact. 
Therefore, we calculate dynamical masses based on both 
a series of spherical models and a disk model. These 
models are also described by Shier (1995). 

Since these systems may not be in dynamical equilib- 
rium, one could question the use of this type of model. 
Clearly, LIRGs with tidal streams and tails are not in 
equilibrium in their outer parts. However, our models are 
generated exclusively for the inner 1 kpc of each galaxy, 
where more stability might be expected and where the 
dynamics have been shown to adhere well to the funda- 
mental plane. Very high resolution molecular line (GO) 
measurements of LIRGs could be used confirm the rela- 
tive equilibrium of the inner parts of the LIRG systems 
and confirm the values of the velocity dispersion obtained 
from the 2.3 /im GO bandheads. 

3.3.1. Spherical Models 

The mass density distributions in the spherical models 
are taken from Tremaine et al. (1994), who analytically 
derive a family of "yy-models" , where the central density 
cusps are described by p oc r^"**, where p is the mass 
density. All 77 models have p oc r^^ as r -^ 00. The 



general form of the density distribution function is 
V 1 
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Models with 77 = 1 resemble a singular isothermal sphere, 
?7 = 2 models have p ~ 1/r at small radii, and rj = 3 
models have constant density cores. We consider only 
the r) = 2 and the rj = 3 models for the rest of this work, 
following Shier (1995). The mass interior to the radius 
is given by 
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We combine these models with the hydrostatic equilib- 
rium equation. 
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where aios is the line-of-sight velocity dispersion and $ 
is the gravitational potential, and solve for the integral 
form of the velocity dispersion. 



(^fosir) = 
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M{r)p{r) dr 



(4) 



However, stars at all different radii in the galaxy con- 
tribute to the observed velocity dispersion within the slit 
area of the spectrograph, in essence taking an average of 
all of the velocities. This average is weighted by the flux 
from each volume element. In terms of the mass-to- light 
ratio, 7, and the extinction along the line-of-sight, e""^, 
the observed velocity dispersion is given in terms of the 
line-of-sight dispersion by 



L<s{r)phe-^dV 
J^p/je-dV 



(5) 



This is a volume integral that covers the portion of the 
LIRG seen through the aperture. 

We can rewrite the form of p in terms of an effective or 
scale radius, Reff, that describes the physical size of the 
model galaxy and the radius at which the density power 
law changes. This form is 



p{r) 
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{r/R,ffyi-^l + r/R,ffy+^ ' 
We can then write the mass interior to a radius r as 
R^ffal {r/R^ff^ Lp/ve'^dV 
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In this case w^ „ (r) is a dimensionless velocity dispersion 
calculated by Tremaine et al. (1994). The definition is 
given by 



<,(r) = GMr/RefMr), 



(8) 



where Mt represents the total mass of the galaxy. Ta- 
ble 3 lists the calculated LIRG masses for 77 = 2 and 
?7 = 3 spherical models based on our measured velocity 
dispersions and scale-lengths. 

3.3.2. Disk Models 

To create a simple disk model, we assume that the local 
velocity dispersion is small in comparison to the circular 



velocity, Vc, so that the stars in the disk are on nearly 
circular orbits, given by 



2_ _ GM{r) 



vi{r) 



(9) 



The inclination reduces the line-of-sight velocity by cosz. 
Also, stars that are not on the observed major axis have 
components of their motion in the plane of the sky, so 
that the observed velocity is actually given by 



Vie 



5 (r, 9) ~ Vc cos i sin 9 



(10) 



where 9 is the azimuthal coordinate in the galaxy. As we 
have no direct information about the inclination angles 
for the LIRGs, we set i equal to 30°, which is the me- 
dian inclination for a set of randomly oriented disks. The 
observed velocity dispersion is the standard deviation of 
the velocities of the stars within one spectrograph aper- 
ture width. This standard deviation is weighted by the 
flux received at each point along the aperture, so that 
the observed velocity dispersion, Ca, is given by 
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(11) 



Assuming that the LIRG has an exponential surface mass 
density, i.e., S = Y^^e^/^'^p, the mass of the galaxy within 
a specific radius is 
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4. COMPARING MOLECULAR MASSES TO DYNAMICAL 
MASSES 

Table 3 shows the molecular gas masses from millime- 
ter CO measurements assuming the standard ^^CO-to- 
H2 conversion factor (Scoville et al. 2000), dynamical 
masses from other works (Shier et al. 1996; Alonso- 
Herrero et al. 2001; Tacconi et al. 2002), and the dy- 
namical masses from this work for the disk, rj = 2, and 
77 = 3 models. All dynamical masses are the values for 
the galaxy mass within the central 1 kpc, for direct com- 
parison with the Shier et al. (1996) values. The mass 
errors are calculated from the errors in the velocity dis- 
persions derived from the near-IR spectroscopy. 

It is interesting to compare our dynamical mass for 
NGC 1614 with that from Alonso-Herrero et al. (2001). 
They used Br 7 spectroscopy (Puxley & Brand 1999) and 
HST/NICMOS images in Pa a of an emitting ring region 
to calculate a dynamical mass independent of the CO 
bandheads. They assumed that the ring is associated 
with a Lindbland resonance or similar dynamical feature 
and that it is close to round. They obtained a rotational 
velocity for the ring by modeling the line profile expected 
for a rotating ring and also calculated an inclination for 
the ring based on its observed ellipticity (i = 51°). They 
then used the virial theorem to convert the ring's rota- 
tional velocity to a mass. They found a dynamical mass 
of 1.3 X 10^ M© within '^620pc. They also estimated 
the mass in old stars as 1.2 x 10^ Mq, and the modeled 
starburst mass as 0.55 x 10^ Mq, both of which fit within 
our mass budget for the stellar population in that region 



in addition to the Scoville et al. (2000) gas mas estimate. 
They noted that their mass determination was lower than 
that of Shier et al. (1996), and this discrepancy is in- 
creased relative to our new determination. Changes in 
the estimate of the inclination angle of the system do not 
increase their dynamical mass by a large enough factor 
to make up the difference in estimates. They use the 
Shier et al. (1996) scale-lengths for their mass models, 
but even adopting those scale-lengths with our value for 
the velocity dispersion does not bring the estimates into 
agreement. This comparison emphasizes the sources of 
error in the dynamical mass determinations. It suggests 
that, at least for some of our sample, our measurements 
should be interpreted as upper limits because the dis- 
persions are measured only on the nuclei and a compact 
nuclear star cluster with very strong CO bands may bias 
our dispersions toward high values. Conversely, a galaxy 
nucleus dominated in K-band by an active nucleus would 
have a bias in dispersion measurements toward anoma- 
lously small values; however, except for Mkn231, none 
of our sample have strong AGN. 

In the cases of NGC 1614 and NGC 2623, the higher 
dynamical mass estimates now allow for the molecular 
gas portion calculated by Scoville et al. (2000). The 
most favorable (largest dynamical mass estimates) cases 
leave ~ 3x 10^ M© for the stellar population in NGC 1614 
and ~ 7 x 10^ M© for NGC 2623. However, for half (5) 
of galaxies in the sample, we find total dynamical masses 
that are lower than the molecular gas masses predicted 
from millimeter CO measurements using the standard 
conversion between CO and H2. The gas masses range 
from ^^1.5 to 10 times the dynamical masses in these 
LIRGs. This is not reasonable, assuming that at least 
some of the mass in the inner regions of these galaxies is 
in the form of stars. For two other galaxies, NGC 3690 
and Mkn 273, almost the entire dynamical mass budgets 
would have to be taken up by molecular gas for the esti- 
mates to agree. 

The results for these ten LIRGs and ULIRGs have 
strong implications for theories on the efficiency and 
mechanisms for star formation in extreme environments 
such as the inner portions of LIRGs. The Schmidt law 
(Schmidt 1959) is widely applied to fit star forming be- 
havior, where the star formation rate is given by a gas 
density power law, 

SsFfl = A^gas 1 (13) 

where A is the absolute star formation rate efficiency. 
Observations generally place the value of N in the range 
1-2, depending on which star formation tracers are used. 
Kcnnicutt (1998) explored the changes in A and N in 
the high density environments of a sample of luminous 
infrared starburst galaxies, including several LIRGs in 
our sample. He found a well-defined Schmidt law for 
the starbursts, with N = 1.3 — 1.4, and calculated that 
the global star formation efficiencies are larger than for 
a sample of normal spirals by a median factor of 6. How- 
ever, he cautioned that the CO/H2 standard conversion 
factor led to scatter in correlations between the star for- 
mation rate and molecular gas density. For instance, it 
may be that the conversion factor is valid in regions with 
solar metallicity but that it underestimates the H2 mass 
in metal-poor galaxies. He found that metal-rich spi- 
rals showed a better defined star formation rate versus 
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molecular hydrogen gas density, in support of the fact 
that the CO/H2 conversion might indeed be wrong for 
some environments. He further showed that lowering gas 
masses in the starburst galaxies by a factor of two, cor- 
responding to a conversion factor half of the standard 
value, increased N from 1.4 to 1.5. With the indications 
that the overestimate of the conversion factor may be 
more like a factor of 3 — 4 than two (this work; Maloncy 
& Black 1988; Solomon et al. 1997; Downes & Solomon 
1998; Lisenfeld, Isaak, & Hills 2000), it seems likely that 
N lies between 1.5 and 2. As a result, the star formation 
efficiency must increase substantially with increasingly 
luminous star formation. 

Our results may also indicate what mechanisms con- 
trol the star formation in merging starbursts. There are 
three basic physical processes behind intense starbursts 
such as those in LIRG environments: gravitational in- 
stabilities in nuclear gas disks (Kennicutt 1998), efhcient 
cloud-cloud collisions driven by gas motion in their cen- 
tral regions which are proportional to a power of local 
volumic density, oc p^ (Scoville et al. 1991), and dy- 
namical disturbances or tidal forces, perhaps, for exam- 
ple, from supermassive binary black holes (Taniguchi & 
Wada 1996). Schmidt laws with N ^ 1 — 1.5 favor the 
gravitational instability scenario, while laws with N '^ 2, 
such as are suggested for our LIRGs, favor star formation 
triggered by cloud-cloud collisions. The third mechanism 
of dynamical disturbances also relies on gravitational in- 
stability formalism, so is generally associated with lower 
values of N. 

Taniguchi & Ohyama (1998) found that, for two sam- 
ples of high luminosity starburst galaxies with millimeter 
CO measurements, the surface density of the far-infrared 
luminosity was proportional to the derived H2 surface 
mass density, with A'^ ~ 1 , implying that the most likely 
star formation mechanism is gravitational instability in 
nuclear gas disks when the standard ratio is used. How- 
ever, Scoville, Sanders, & Clemens (1986) suggested that 
cloud-cloud collisions account for the high rates of star 
formation in LIRGs, with this process being the domi- 
nant mode for the formation of high-mass stars. They 
explored CO data on a large sample of giant molecu- 
lar clouds and suggested that any model attempting to 
mimick massive star formation must account for both 
the formation of stars deep in the interior of clouds and 
the lower efficiency in high-mass clouds. A model that 
accounted for both of these factors was one in which the 
star formation is triggered by cloud-cloud collisions, so 
that OB star clusters appear on the edge of and near the 
center of merged clouds. The lower gas masses implied 
by the calculated dynamical masses for our sample sug- 
gest larger values of A^ , which can only be supported by 
the more efficient, high mass star formation allowed by 
cloud-cloud collisions. 

Based on our dynamical masses and the works cited 
above, we conclude that at least some LIRGs do not 
harbor the enormous molecular gas masses calculated 
from CO emission measurements and that the standard 
conversion factor from CO to H2 derived from local gi- 
ant molecular clouds cannot be used in the dense in- 
ner environments of starburst galaxies. Our dynamical 
masses essentially place upper limits on the amount of 
molecular gas likely to be at the centers of LIRGs. We 
further predict that the star formation efficiencies are 



higher than those derived from the standard conversion 
values and that the dominant source of star formation in 
LIRGs is likely to be cloud-cloud collisions. We expect 
that LIRGs do not have prohibitively high extinctions 
at near-infrared wavelengths and anticipate that much 
of their structure and starburst information is well diag- 
nosed with H and AT^-band observations. This conclu- 
sion is supported by mid-infrared observations such as 
Gallais et al. (2004) where Ay is found to be on the 
order of tens rather than hundreds. 

5. SUMMARY 

We have explored the dynamics of a sample of ten lu- 
minous and ultraluminous infrared galaxies using near- 
infrared imaging and spectroscopy. Velocity dispersions 
for the galaxies were calculated from the CO absorption 
bandheads present at 2.3/im. Exponential disk and effec- 
tive bulge scale-lengths for each LIRG were derived using 
two-dimensional model fits to i/-band images. These re- 
sults were used in conjunction with a series of spherical 
and disk-like density profile models to determine total 
dynamical masses for the galaxies. 

The dynamical masses for the LIRGs are in the range 
^ 2 X 10^ — 4 X 10^° M© for the innermost 1 kpc. We com- 
pare these masses to estimates of molecular gas mass de- 
rived from interferometric millimeter wavelength obser- 
vations of CO (e.g., Bryant & Scoville 1999). The molec- 
ular gas masses are usually calculated assuming a stan- 
dard conversion between CO luminosity and the amount 
of H2, based on ratios found for giant molecular clouds 
in the Milky Way. This conversion has come under some 
criticism from both theoretical (Maloney & Black 1988) 
and observational (Shier, Rieke, & Rieke 1996) points 
of view, and it has been suggested that using such a 
conversion in the extreme environments of the nuclei of 
starburst galaxies overestimates the amount of molecular 
hydrogen present. In agreement with these criticisms, we 
find that for over half of the LIRGs in our sample, the 
molecular gas masses exceed or fill the entire dynami- 
cal mass budget, an unreasonable result considering the 
stellar luminosity that is observed. We conclude that, 
at least for some LIRGs, the standard conversion is an 
inappropriate tool for deriving molecular gas masses. 

These results have wide implications for the efficiency 
of star formation in LIRGs and for the physical mech- 
anisms responsible for such star formation. Decreased 
amounts of molecular gas, as implied by our dynami- 
cal masses, allow an increase in the power of the global 
Schmidt law, suggesting that the efficiency and rate of 
star formation in LIRGs is higher than predicted by mil- 
limeter observations alone (Kennicutt 1998). Such large 
star formation rates are thought to have been created by 
efficient cloud-cloud collisions in the inner regions of the 
LIRGs, dependent on the square of the density in the 
area (Scoville, Sanders, & Clemens 1986), rather than 
by gravitational instabilities in nuclear disks or other 
dynamical or tidal disturbances (Taniguchi & Ohyama 
1998). We also conclude that the huge amounts of ex- 
tinction implied by large molecular gas masses may not 
be present and that near-infrared observations are often 
adequate probes of the dynamics and starburst activity 
associated with LIRGs. 
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TABLE 1 
The LIRG Sample and Observations 



Name 
(1) 


Morphological 
Type 

(2) 


(Lq) 

(3) 


Imaging 
Band 

(4) 


Exposure 

Time (min) 

(5) 


Imaging 
Date 

(6) 


Instrument/ 
Telescope 

(7) 


Spectroscopy 
Date 

(8) 


Exposure 

Time (min) 

(9) 


NGC 1614 


B(s)c pec; Hll: Sy2 


3.8x10" ^ 


F160W 


3 


1998 Feb 7 


HST/NICMOS 


2000 Dec 5 


60 


NGC 2623 


merger 


3.5x10^^ ^ 


H 


21 
19 


2000 Dee 9 
2000 Dee 9 


PISCES/2.3 m 


2000 Dec 5 


108 


NGC 7674 


S(r)bc pee; Hll: Sy2 


3.1x10" ^ 


H 


21 
20 


2000 Dee 9 
2000 Dee 9 


PISCES/2.3 m 


2000 Dec 5 


48 


IRAS 10565+2448 


Hll 


10^^ a 


H 


13 

8 


2001 May 8 
2001 May 8 


PISCES/2.3 m 


2000 Dec 5 

2001 Mar 2 


36 

72 


IRAS 17208-0014 


Sbrst H II 


2.5x10^^ "^ 


H 


10 
10 


2003 Mar 13 
2003 Mar 13 


PISCES/6.5 m 


2001 Mar 2 
2001 Mar 3 


12 
36 


IC 694/NGC3690 


Sbrst AGN 


5.2x10" " 


H 


10 
10 


2001 Apr 9 
2001 Apr 9 


PISCES/6.5 m 


2001 Mar 3 


84 


IC 883 


Im:pee; Hll LINER 


4.0x10" <= 


F160W 


6 


1997 Nov 21 


HST/NICMOS 


2001 Mar 3 


60 


VII Zw31 


Hll 


8.7x10" '^ 


H 


7 


2001 May 8 


PISCES/2.3 m 


2001 Dec 5 


144 


Mkn 231 


SA{rs)e? pec; Syl 


3.5x10^^ <= 


H 


10 

7 


2001 May 8 
2001 May 8 


PISCES/2.3 m 


2001 Mar 2 


48 


Mkn 273 N 


Sy2 


1.4x10^^ •= 


H 


10 


2003 Mar 11 


PISCES/6.5 m 


2001 Mar 3 


84 



Note. — Column header explanations. Col. (1): Gala 
(3); IR luminosity (Lq) from 8-1000 ^m. Col. (4): Nea 
Col. (7): Instrument and telescope used for the near-inf 



e. Col. (2): Morphological type from the NASA Extragalactic Database or Scoville ct al. (2000). Col. 
cd imaging band. Col. (5): Imaging exposure time in minutes. Col. (6): Date of imaging observation. 
Tiaging. Col. (8): Spectroscopy date of observation. Col. (9): Spectroscopy exposure time in minutes. 



L et al. (1995) 



Ki 

Gallais et al. (2004) 

'^Sanders et al. (1988) 

Scoville et al. (1989) 



TABLE 2 
LIRG Scale-lengths and Velocity Dispersions 



Name 


Rexp 


-Re// 


CF 


(Mother 




(pc) 


(pc) 


(kms-1) 


(kms-1) 


(1) 


(2) 


(3) 


(4) 


(5) 


NGC 1614 


85.8 ±58.9 


58.9 ±105.1 


164 ±8 


75 ± 12=^ 


NGC 2623 


1883.0 ±399.8 


1007.8 ± 149.2 


153 ±11 


95 ± 13=^ 


NGC 7674 


4420.5 ±215.0 


890.8 ±181.5 


93 ±26 




IRAS 10565+2448 


1454.3 ± 643.6 


455.5 ±87.8 


141 ±4 


125 ±31*^ 


IRAS 17208-0014 


1618.6 ±25.9 


932.1 ±73.2 


125 ± 28 




IC 694/NGC 3690 


686.5 ±28.9 


12.15 ±4.19 


141 ±17 


135 ±21=^ 


IC 883 


224.6 ± 18.4 


19.0 ±48.9 


151 ±5 




VII Zw 31 


1533.3 ±53.4 


416.7 ±320.6 


98 ±13 




Mkn 231 


573.6 ±21.1 


400.7 ±46.4 




115 ±10^= 


Mkii 273 N 


1317.0 ±10.8 


732.4 ±82.0 


232 ± 43 


285 ± 30^= 



Note. — Column header explanations. Col. (1): Galaxy name. Col. (2): iJ-banc 
tial disk scale-length in parsecs. Col. (3): ff-band effective radius of the bulge i 
Col. (4): Velocity dispersion in km s~^ from this work. Col. (5): Velocity dispen 
s~lfrom other works. 

^Shicr ct al. (1996) 

Dasyra ct al. (2006) 
■^Tacconi ct al. (2002) 
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TABLE 3 










LIRG Masses 






Name 


logM(H2) 


log Mdyn 


(Mq) 


(Mq) 


(Mq) 




(Mq) 


(Mq) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC 1614 


9.68 


9.32+0.25=^ 


10.01±0.05 


9.72±0.06 


9.84±0.07 


NGC 2623 


9.77 


9.46±0.2l'' 


10.00±0.06 


10.22±0.07 


10.11±0.07 


NGC 7674 


10.12 




9.52±0.29 


9.82±0.28 


9.73±0.28 


IRAS 10565+2448 


10.34 




10.03±0.10 


10.16±0.10 


10.17±0.10 


IRAS 17208-0014 


10.71 




9.92±0.22 


10.09±0.22 


9.99±0.22 


IC 694/NGC 3690 


10.00 


9.76±0.24b 


lO.OliO.ll 


10.13±0.11 


10.09±0.11 


IC 883 


9.87 




10.09±0.03 


9.22±0.03 


9.29±0.03 


VII Zw 31 


10.69 




9.73±0.12 


9.60±0.12 


9.09±0.12 


Mkii 231 


10.33 


10.60^= 


9.90±0.08 


9.97±0.08 


10.00±0.08 


Mkn 273N 


10.58 




10.47±0.18 


10.62±0.18 


10.57±0.18 



Note. Column header explanations. Col. (1): Galaxy name. Col. (2): Log of the molecular hydrogen 

mass from ScoviUe et al. (2000) in solar masses, assuming the standard CO-to-H2 conversion factor. Col. 
(3): Log of the dynamical mass in solar masses from other works. Col. (4): Log of the calculated dynamical 
mass in solar masses from this work for a disk model. Col. (5): Log of the calculated dynamical mass in 
solar masses from this work for a spherical, rj = 2 model. Col. (6): Log of the calculated dynamical mass 
is solar masses from this work for a spherical, rj = 3 model. 



Alonso-Hcrrcro ct al. (2001) 
°Shicr ct al. (1996) 
■^Tacconi ct al. (2002) 
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NGC1614 



Mac 1621 



NGC7674 



---ir- -" •^J: 








Fig. 1. — The _ff-band images of the LIRG sample with a field of view of approximately 3' X 3', with north up and east to the left. The 
two HST/NICMOS images (NGC 1614 and IC883) have fields of view of ~30" X 30". 
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Fig. 2. — Onc-dirncnsional spectra of the LIRG sample showing the CO absorption bandheads caused by the 2-0 and 3-1 molecular 
transitions. Mkn231 was observed at 1.67 ^m to measure the second overtone CO bands because of the dilution of the first overtone bands. 



